Domain-wall complexes in ID ferromagnets and critical media 
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Interactions of domain walls (DWs) are analyzed with relevance to formation of stationary bubbles 
(complexes of two DWs) and complexes of many domains in one dimensional systems. I investigate 
the domain structures in ferromagnets which are described with the Landau-Lifshitz equation as well 
as the domains in critical systems described with the Ginzburg-Landau equation. Supplementing 
previous author studies on the creation of hard bubbles [formed by one Bloch DW and one Neel 
(Ising) DW] in the presence of an external (magnetic) field, the soft bubbles consisting of two Bloch 
DWs or two Neel (Ising) DWs are studied in detail. The interactions of two DWs of the same kind 
are studied in the framework of a perturbation calculus. 
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I. INTRODUCTION 

Localized and patterned (labirynth) structures in 
bistable media are widely considered with relevance to 
the storage of binary information [l|, Q ■ Far enough from 
the phase-transition point, the complexes of DWs are 
observed in chemical reactors, magnetic and polar (solid 
and liquid) media while their basic properties are de- 
scribed with equations of motion with are similar for dif- 
ferent media. In particular, far from the critical regime, 
the magnetization (polarization) dynamics is described 
with the Landau-Lifshitz-Gilbert (LLG) equation while 
near the criticality it is governed by the Ginbzurg-Landau 
(GL) equation. Recently, one observes an especial inter- 
est in ID complexes of magnetic and polar DWs due to 
hopes for miniaturization of memory carriers which result 
from achievements of current nanowire technology |3- 5] . 
The problem of stability of many-domain complexes is 
connected to the need of switching the domain-encoded 
binary information on demand since such process induces 
unbalanced interactions of the DWs. 

In the present work, I analyze binary interactions be- 
tween the DWs (of the Bloch and Neel-Ising type) with 
dependence on the distance of their separation and their 
chiralities (opposite or like) and I study the formation of 
magnetization bubbles in ID far from and near the crit- 
icality. It enables me to consider the stability of many- 
domain structures. With correspondence to 2D magnetic 
bubbles (which are called hard ones when their bound- 
ary is composed of alternating Bloch and Neel lines), I 
divide the ID bubbles into hard ones (composed of one 
Bloch and one Neel-Ising DWs) and soft ones (composed 
of two Bloch DWs or two Neel-Ising DWs) @,0|< In my 
recent papers on the externally-driven motion of DWs, I 
have studied the formation of ID hard bubbles via colli- 
sion of a Bloch DW with a Neel (Ising) DW in the normal 
and critical regimes, (solving the LLG and GL equations, 



respectively) [1, @. Here, the energy of soft bubbles is 
studied as function of the bubble length within a per- 
turbation calcus in order to complete the picture of the 
DW interactions. I follow a perturbation approach to 
the DW interaction developed in Ref. [10( with relevance 
to the parametrically-driven nonlinear Schrodinger equa- 
tion. It differs from previous ones (e.g. [TTj] ) in terms of 
the perturbation expansion of the dynamical parameter 
(magnetization) whose present form ensures conservation 
of the magnetization length. 

With application to ferromagnets, the ID approxi- 
mation is relevant to crystalline nanowires with strong 
bulk anisotropy compared to surface magnetostatic ef- 
fects (and, especially, with circular cross section). The 
present analysis provides a basis for comparison of the 
DW interactions in such magnets to the interactions of 
DWs in the nanostripes of noncrystalline ferromagnets 
which are studied in Ref. 

Binary interactions od the DWs and the stability of 
ID bubbles are considered in sections II and III, with 
relevance to the ferromagnetic wire at low temperatures 
and to the subcritical systems, respectively. In section 
IV, the stability of many-DW ID systems is considered 
on the basis of previous-section results. 



II. MAGNETIC BUBBLES AT LOW 
TEMPERATURES 

I consider DW solutions to the LLG equation 
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The first term of the r.h.s. of (fTJ) relates to the exchange 
spin interaction while the second (Zeeman) term depends 
on the external magnetic field H = (H x ,0,0), thus, 7 
denotes the giromagnetic factor. The constant Pi(2) de- 
termines strength of the easy axis (plane) anisotropy and 
i = (1, 0, 0), j = (0, 1, 0). I notice that the long axis of a 
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nanowire is an easy axis for the most often investigated 
magnets, however, another choice of the anisotropy axes 
does not influence the magnetization dynamics. Since 
LLG equation is valid only when the constraint |m| = M 
is satisfied, one writes equivalent to ([1]) equations of the 
unconstrained dynamics. Introducing m± = m y ± im z , 
one represents the magnetization components using a 
pair of complex functions g(x,t), f(x,t). The relation 
between the primary and secondary dynamical variables 
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ensures that |m| = M. Inserting ([2]) into ([TJ leads, fol- 
lowing the Hirota method for solving nonlinear equations 

HQ, to 
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where A, A denote Hirota operators of differentiation 
which are defined by 

D\ n D n x b(x, t) ■ c(x, t) = (d/dt - d/dt) m 
x [8/dx - d/dx') n b(x, t)c(x , i)\ x=x >, t=t > . 

For H = 0, the stationary single-DW solutions to (|3]) 
take the form 
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where \k\ = y/ j3i / J, w = —w*, (a Bloch DW) or |fc| = 
/3 2 )/J, w = w*, (a Neel DW). Let us define ip 
following e lip = w/\w\, thus, <p — ±ir/2 for Bloch DW 
while <p = 0, 7r for Neel DW. 

Since neither exact double-Bloch nor double-Neel so- 
lutions to the LLG equation are not known for the case 
of zero external field, (the length of the soft bubbles di- 
verges with H x — > 0, jla - fr7| ). I analyze the interactions 
of the relevant DW pairs (the pair of Bloch DWs and the 
pair of Neel DWs) within a perturbation calculus. Lo- 
cally, in the vicinity of the center of jth DW, (j — 1,2), 
one can write the magnetization in the form 



m(ir, 0) = m 
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where m^' denotes a stationary single-DW solution to 
(p} [which corresponds to (J4j) ] 



m+\x) = Me iVj sech[(Tjk(x — xoj)], 

m xH x ) — —Mtaiih[ajk(x — xoj)] (6) 

while Smy) denotes a perturbation due the presence of 
another DW. Here o\ — —<t 2 — ±1, ip\ — + f 2 and \y>j \ = 
7r/2 (a pair of Bloch DWs) or tpj — 0, tt (a pair of Neel 



DWs). When assume k > 0, Uj = 1 relates to the head- 
to- head spin structure while <Zj = —1 to the tail-to-tail 
structure. In order to satisfy the constraint |m| = M, I 
apply the perturbation of the form 
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Plus and minus relate to the bubble magnetization par- 
allel and antiparallel to x-axis, respectively. This form 
of the initial magnetization reflects the fact that the in- 
teraction of topological solitons (DWs) is not simply due 
to their overlap while it is accompanied by some reorien- 
tation of whole the separating them domain. Therefore, 
the magnetization © contains only products of the com- 
ponents of m' 1 ) and m^ 2 ) while it does not contain linear 
in m^) terms. Let a = k(xo 2 — xqi). Inserting (|SJ) into 
the Hamiltonian 
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where Hz denotes the Zeeman part of the Hamiltonian, 
one arrives at the dependence of the energy E (z) = 
Ho(z)dx on the distance between the DW centers 
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where 9 = 1 for Neel DWs while 9 = for Bloch DWs, 
It 



and Ir (a) denote integrals given in Appendix. The upper 



signs correspond to the pair of the DWs of the opposite 
chiralities, the case ipi = <f2, while the lower signs cor- 
respond to the pair of the DWs of like chiralities, the 
case cp i = if2 + tt- According to the plot of energy of 
the DW pair (Fig. 1), in absence of the external field, 
the interaction is attractive when both the DWs are of 
opposite chiralities while it is repulsive in the case of like 
chiralities. 

In both the cases (of opposite and like chiralities), it 
is possible to create a stationary magnetization bubble 
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FIG. 1. Energy of a pair of DWs as a function of the distance 
of their separation; solid (dash-dotted) line - DWs of opposite 
chiralities (an untwisted pair) in absence (presence) of an ex- 
ternal field, dashed (dotted) line - DWs of the same chiralities 
(a twisted pair) in absence (presence) of an external field. 



applying an external magnetic field in the direction pur- 
suant or opposite to the bubble magnetization, respec- 
tively. However, only in the case of ifi = ip% ± 7r, the 
function E(a) = Eq(cl) + Ez(a) has a minimum at a =/= 0, 
(see Fig. 1). Therefore, only the bubble created by the 
DWs of like chiralities is stable against external pertur- 
bations. This stable bubble corresponds to the final state 
of the long-term evolution of a breather. 

For H 7^ 0, one can find strict double-Bloch and 
double-Neel solutions to LLG equation, thus, one can 
verify our prediction on existence of soft stationary bub- 
bles. Assuming the stationary double-Bloch DW solution 
and double-Neel DW solution to ([3]) to be of the form 
/ = 1 + v. g = we kx + we~ kx , v = v* , one finds 
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for w = to*, 9 = 1 (an untwisted double-Neel DW), and 
for w — — w* , 9 = (an untwisted double-Bloch DW). 
Inserting / = 1 +v, g = we ' 



we into © , one finds 
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for w — w* , 9 — 1 (a twisted double-Neel DW), and 
for w = — w*, 9 = (a twisted double-Bloch DW). 
The untwisted double-DWs are called nuclei while the 
twisted double-DWs are called 27T-DWs. Since typically 
Pi S> P2 while \H X \ = /3i J '7 corresponds to the coercion- 
field value, (hence, I assume \H X \ <C Pi/j), the untwisted 
and twisted double-DWs relate to different (parallel or 
antiparallel to the bubble magnetization, respectively) 
directions of the external magnetic field, which ensures 
that v in (1121) or (ITBl is determined. In order that the 




FIG. 2. Formation of hard bubble (composed of one Neel DW 
and one Bloch DW) in a longitudinal field. The DW reflection 
is accompanied by a change of the Bloch wall (of the velocity 
V2) into the Neel wall (of the velocity v\) and vice versa. 



existence of nuclei was possible, the DWs of like chirali- 
ties must attract, thus, the nucleus vanishes after turning 
the magnetic filed off. The prediction on stability of the 
27r-DWs and instability of nuclei is in accordance with a 
study of linear-excitation scattering on the soft magneti- 
zation bubbles and with simulations [TBI . [l7j . 

In order to complete the picture of the bubble forma- 
tion in ID ferromagnet, I refer to my previous study of 
complexes of one Bloch DW and of one Neel DW Q. In 
the absence of external field, such a stationary complex is 
described by a strict solution to the LLG equation, hence, 
there is no interaction between the walls. An interaction 
appears, however, upon the field application due to a 
dynamically induced deformation of the DWs. The ex- 
ternal field enforces motion of any DW and its anti-wall 
in the opposite directions and it can induce the collision 
of a Neel DW with a Bloch DWs. The result of such a 
collision has been found to be their mutual reflection ac- 
companied by the transformation of the Bloch wall into 
the Neel DW and vice versa. The process can be consid- 
ered as the transmission of the Bloch DW through the 
Neel DW with change of the head-to-head structure into 
the tail-to-tail one (and vice versa) which is illustrated in 
Fig. 2, and it is similar to the collision of spontaneously 
propagating (in absence of dissipation) DWs |l4|. Since, 
under the action of the external field, the different-type 
DWs interact repulsively, they can form a bubble which is 
an ID counterpart of the hard bubble in planar magnets. 



III. BUBBLES NEAR CRITIC ALITY 

The subcritical dynamics of bistable systems is gov- 
erned by the GL equation 
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For Pi > 3/?2, (HH describes the evolution of Bloch (of 
lower energy) and Ising (of higher energy) DWs and their 
complexes. The field-induced creation of ID hard bub- 
ble (composed of one Bloch DW and one Ising DW) has 
been investigated in Q. It has been found to be similar 
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to the dynamics of the hard bubble in the low tempera- 
ture regime. For H — 0, the stationary complex of one 
Ising DW and one Bloch DW in an ID subcritical system 
is described by a strict solution to (|14p . thus, both the 
DWs do not interact [llj]. The field-driven collision of 
the Ising DW with the Bloch DW induces their mutual 
reflection, hence, they can form a similar hard bubble 
as it was described in the previous section. I mention 
that an analysis of linear-wave scattering on the Bloch- 
Ising complex with the parametrically-driven nonlinear 
Schrodinger equation has shown the dynamically-induced 
repulsion of both the walls also [19J ■ Below I pay my at- 
tention to soft magnetization bubbles (of two Bloch DWs 
or of two Ising DWs). 

In the first perturbation approximation, in the vicinity 
of the center of jth DW (j=l,2), the magnetization profile 
can be written as 



m(x, 0) = m^' (x) + SmP' (x), 



(15) 



where denotes a strict (stationary) single-DW solu- 
tion to (Hill for H = 



-»w - ^ 
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Here a\ = — a 2 = ±1, and k = \/2j3 2 / J, <Pj — ±n/2 
(Bloch DWs) or k = y/ljh + /3 2 )/(2 J), ipj = 0, vr (Ising 
DWs). According to my claim on the form of the 
ansatz (j8|), since the interaction of DWs is related to 
a reorientation of the separating them domain, I ex- 
pect the perturbed magnetization not to contain lin- 
ear in mP' terms. Taking the following perturbation 
5 m W = fflW[W^/i/(ft +7h)- 1], (fc ^ j), leads 
to 
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This form of the initial magnetization was used in [To| . 
whereas, previous treatments of the DW interactions 
within the GL approximation used the ansatz Sm^ — 
± yj {Pi + f3 2 )/fi, (k 7^ j), by an analogy to the 
perturbation calculus for nontopological (bright) solitons 

Inserting (JTTj) into the Hamiltonian 
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one arrives at the dependence of the energy E (z) = 
J^o^o(z)da; on the distance between the DW centers. 
Defining a = k(xo 2 — xqi), for the interacting Ising DWs, 
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FIG. 3. Energy of a pair of DWs as a function of the distance 
of their separation; a) solid (dash-dotted) line - Ising DWs in 
absence (presence) of an external field, b) solid line - Bloch 
DWs of the opposite chiralities (an untwisted pair), dashed 
(dotted) line - Bloch DWs of the same chiralities (a twisted 
pair) in absence (presence) of an external field. 



one finds 
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while for the Bloch DWs, (Eq is derived up to the lowest 
order in /3 2 //3i); 
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The integrals Ifc(a) are given in Appendix. Here, the 
upper signs correspond to a pair of the Bloch walls of 
opposite chiralities, the case ip\ — ip 2 while the lower 
signs to a pair of the Bloch walls of like chiralities, the 
case ipi = if2 ± vr. Plotting the energy Eo(a) of the DW 
pairs (Fig. 3), we see the Ising DWs to attract each other 
in the absence of external field while the character of the 
interaction of Bloch DWs to be dependent on their chi- 
ralities. In the presence of an external field parallel to 
the bubble magnetization, the Ising walls can form a sta- 
tionary state corresponding to the maximum of energy, 
however, it is ustable to perturbations. The Bloch DWs 
of opposite chiralities attract each other when the dis- 
tance of their separation is big and they repel each other 
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at short separation distances. Because of the minims of 
the energy E (a) at a ^ 0, stable bubble of Bloch DWs 
(a breather) can be formed even in absence of the field. 
The Bloch DWs of like chiralities always repel in absence 
of the field, hence, an external field antiparallel to the 
magnetization of the intermediate domain enables cre- 
ation of a stable bubble (a breather) which corresponds 
to the minims of Eq(cl) + Ez(a). 

IV. STRUCTURES OF MANY DWS 

Stability of many-DW structures is an important prob- 
lem in view of the challenge of designing a nanowire-based 
information carrier (a DW-racetrack memory) [3[. Peri- 
odic structures of the Bloch (Neel) DWs are stationary 
because of the compensation of (attractive or repulsive) 
interactions between the walls. However, bit recording 
requires switching the magnetization of a single domain 
of the memory tape on demand. For instance, the magne- 
tization reversal can be performed via flip of the chirality 
of two neighboring DWs in the tape, which initiates their 
movement towards or outwards each other and, eventu- 
ally, results in their annihilation due to the appearance 
of uncompensated attractive forces. Unfortunately, the 
resulting nonperiodicity of the system causes destabiliza- 
tion of the record. I mention that similar process is re- 
sponsible for annihilation of DWs during strong-current 
driven motion of multi-domain systems [201 ] . Then, the 
unbalanced interactions appear because, for the current 



intensity exceeding the Walker-breakdown value, magne- 
tization in different (head-to-head and tail-to-tail) DWs 
rotates in opposite directions about the magnetic axis, 
(unlike upon the application of a strong magnetic field). 

An exception is a subcritical system of Bloch DWs 
whose all the neighboring DWs are of opposite chirali- 
ties and the distance of their separation corresponds to 
the size of the breather found in section III. The result 
of changing the chirality of a single DW or a DW pair 
is some shift of the wall positions without any flip of the 
domains. 

Aperiodic trains of DWs of like chiralities (unstable to 
DW interactions) as well as trains of alternating Bloch 
and Neel (Ising) DWs (unstable to the Neel-Bloch or 
Ising-Bloch transition inside the DWs, [2l|) can be sta- 
bilized by an external magnetic field whose application 
results in creation of soft and hard magnetization bub- 
bles, respectively. The energy of a soft bubble increases 
with the field value [see minims of E(a), dotted lines in 
Fig. 1 and Fig. 3], hence, unlike for weak fields, in a 
regime of strong field, the hard bubbles can be energeti- 
cally favorable [6[. 
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Appendix A: Explicit form of integrals 



( {sech 2 (y)tanh(— y + a) — sech 2 (— y + a)tanh(y) =p [ — tanh(y) + tanh(— y + a)] sech(y)sech(— y + a)} 

-oo 

+ {— sech(y) [tanh(y)tanh(— y + a) + sech 2 (— y + a)] ± sech(— y + a) [tanh(y)tanh(— y + a) + sech 2 (y)] dy 

(Al) 

if (a) = / f[sech(y)tanh(-y + a) ± tanh(y)sech(-y + a)] 2 \ dy (A2) 

J —oo 

pOO 

if (a) = / [tanh(y)tanh(— y + a) =F sech(y)sech(— y + a) + 1] dy (A3) 



h(a) = ( [sech 2 (y)tanh(-y -f a) - sech 2 (-y + a)tanh(y)] 2 + [tanh 2 (y)tanh 2 (-y + a) - l] 2 | dy (A4) 

J —oo 
poo 

I 5 (a)= / [tanh(y)tanh(-y + a) + 1] dy (A5) 
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